Distribution of local Lyapunov exponents in spin-glass dynamics 
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We investigate the statistical properties of local Lyapunov exponents which characterize magnon 
localization in the one-dimensional Heisenberg-Mattis spin glass (HMSG) at zero temperature, by 
means of a connection to a suitable version of the Fokker-Planck (F-P) equation. We consider 
the local Lyapunov exponents (LLE), in particular the case of instantaneous LLE. We establish a 
connection between the transfer-matrix recursion relation for the problem, and an F-P equation 
governing the evolution of the probability distribution of the instantaneous LLE. The closed-form 
(stationary) solutions to the F-P equation are in excellent accord with numerical simulations, for 
both the unmagnetized and magnetized versions of the HMSG. Scaling properties for non-stationary 
conditions are derived from the F-P equation in a special limit (in which diffusive effects tend to 
vanish), and also shown to provide a close description to the corresponding numerical-simulation 
data. 
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I, INTRODUCTION 

The analytic treatment of quenched disordered sys- 
tems in Condensed Matter Physics invokes many con- 
cepts from statistical theory. Among these, we shall be 
concerned in this paper with the connection between the 
probability distribution of Lyapunov characteristic expo- 
nents for the (equilibrium) problem of low-lying magnetic 
excitations in spin glasses, and the Fokker-Planck (F-P) 
equation [3], which is a key element in the description of 
non-equilibrium stochastic processes. 

Links between an F-P equation and stationary distri- 
butions of physical quantities in a transfer-matrix (TM) 
description have been established, usually in the con- 
text of calculating the probability distribution functions 
(PDF) of electronic conductance, in one-dimensional (Id) 
or quasi-ld noninteracting electronic systems with disor- 
der (Anderson localization) @,S,S,H,E0,[i,[il- Specif- 
ically, one considers the distribution of the TM itself (or, 
equivalently, its transmission eigenvalues), between the 
left and right extremes of a wire of length L ^> 1 lat- 
tice spacings. The F-P equation is set up to account for 
the infinitesimal changes in the TM, caused by a small 
length variation SL. The quantities whose PDF is cal- 
culated are, therefore, aggregate in the sense that they 
explicitly incorporate all contributions for the TM, say, 
from x = to L, i.e., they represent global Lyapunov ex- 
ponents. This is in contrast with the approach we take 
here; as shown below, we shall concentrate on the PDF 
of local Lyapunov exponents, which appear not to have 
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been as thoroughly exploited as their global counterparts 
(at least in the Condensed-Matter context of localization 
and similar problems). 

As a model system to apply the ideas developed here, 
we consider spin waves (magnons) in a simplified spin 
glass model, the so-called Heisenberg-Mattis spin glass 
(HMSG) [Tq]. In Id, the dynamics of HMSGs turns out 
to exhibit many non-trivial features fll Jl ^ 13 ^ Til. HH. Ilti 

.[T^|. including dynamic exponents |l5l"lcT l8l d ifferent 
from the standard hydrodynamic predictions (l9l |. 

In Id, all eigenstates are localized for any amount of 
disorder, though as energy uj — > 0, the localization length 
A(w) diverges continuously. This, in turn, suggests the 
applicability of scaling concepts in the low-frequency, 
long- wavelength limit. 

Operationally, one can obtain A(cj) in one-dimensional 
systems by using a TM approach, and extracting the 
smallest Lyapunov exponent (whose inverse corresponds 
to the largest localization length) which arises from re- 
peated iteration of the TM. This has been done for 
HMSG chains (^MEl, and can be extended to HMSG 
in d = 2 and 3~[l3P 

Lyapunov exponents are well defined quantities in the 
sense that, for N 3> 1 iterations of the TM, the width of 
the distribution of their corresponding estimates shrinks 
to zero [23| (as N -1 / 2 in many cases of interest, as a 
consequence of the central limit theorem). However, it 
has been shown 0, [H, SI, S3, S3 that the so-called local 
Lyapunov exponents (LLE, to be defined more accurately 
below) exhibit non-trivial distributions whose width re- 
mains finite even as TV becomes fairly large (compared, 
e.g., with the localization length). Furthermore, LLEs 
may provide a wealth of information specific to the dy- 
namics of the associated physical system, which does not 
manifest itself in the aggregated, <5-function like behavior 
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of their global counterparts. 

In this paper, we investigate the statistical properties 
of local Lyapunov exponents in the Id HMSG at zero 
temperature, by means of a connection to a suitable ver- 
sion of the F-P equation. 

In Section [n] we recall pertinent aspects of the HMSG, 
specializing to Id. In Section Mil we first consider the 
aggregate effects which characterize the (global) Lya- 
punov exponents for the system under study; then we 
investigate local exponents (LLE), in particular the case 
of instantaneous LLE. We establish a connection be- 
tween the TM-iterated recursion relation for the prob- 
lem, and a suitable F-P equation governing the evolution 
of the probability distribution of the instantaneous LLE. 
We show that the closed-form (stationary) solutions to 
the F-P equation are in excellent accord with numerical 
simulations, for both the unmagnetized and magnetized 
versions of the Id HMSG. Scaling properties for non- 
stationary conditions are derived from the F-P equation 
in a special limit (in which diffusive effects tend to van- 
ish), and also shown to provide a close description to 
the corresponding numerical-simulation data. Finally, in 
section [TV} concluding remarks are made. 



II. HEISENBERG-MATTIS SPIN GLASSES: 
SCALING IN ID 

We consider Heisenberg spins on sites of a hypercubic 
lattice, with nearest-neighbor couplings: 



H 



(1) 



The bonds are randomly taken from a quenched, binary 
probability distribution, 

P(J y ) = P 5(J ij - Jo) + (1 -p) S{Ja + Jo) ; (2) 

here we shall mainly consider p = 1/2 (unmagnetized 
spin glass). 

The Mattis model ascribes disorder to sites rather than 
bonds (Jij — > Jo Ci Ci); so that the Hamiltonian reads: 



% — - Jo Ci Cj Si • Sj , 



where Ci = +1, —1 with probability p, (1 — p). Then p = 
1/2 is the Mattis spin glass while p ^ 1/2 corresponds to 
a magnetized Mattis model. This way, the overall energy 
is minimized by making Sf — Q S, which constitutes 
a (classical) ground state of the Hamiltonian Eq. (J3j) . 
Consideration of the spin-wave equations of motion (see, 
e.g., Refs.EHtil) gives, with h = 1: 



iQ dui/dt = V" Jo (ui - uj) 
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(4) 



where the Ui are Mattis-transformed local (on-site) spin- 
wave amplitudes, and the sum is over sites j which are 



nearest neighbors of i. For the eigenmodes with fre- 
quency lu (in units of the exchange constant Jo), Eq. (@]) 
leads to 



LU 



£j = (Ui - Uj 



(5) 



The relationship of frequency to wave number, k (in the 
context of localization, this corresponds to A -1 ), at low 
energies is characterized by the dynamic exponent z: 

lu oc k z . (6) 

For p = 1/2 in Id, it was predicted analyticall y fl5Lll6j |. 
and verified by numerical calculations [20, [2l|, [22| , that 
z = 3/2. Still in Id, but at general p, Eq. ((5j) becomes 



(2 - Ci Ui = Ui-! + u i+1 



(7) 



A TM approach [H Ed] can be formulated, giving 

nana 

MM- 



Ui+l 
Ui 



2-&UJ -1 
1 



Ui 
Ut-l 



Ti{u) 



Ui 



(8) 

The procedure for calculating Lyapunov exponents in 
this case is the same as that used for Anderson local- 
ization problems [2§|. Indeed, in both cases the TM 
is symplectic, and one can use Oseledec's theorem and 
dynamic filtration [23| to extract the smallest Lyapunov 
exponent, whose inverse is the largest localization length. 



III. LYAPUNOV EXPONENTS 

A. Aggregate effects: the TV — > oo limit 

The Lyapunov exponent emerging from iteration of 
Eq. ((HI) is given by 



y[uj) = lim — In 



N 



(9) 



(3) where \v Q ) 



is an arbitrary initial vector of unit 



modulus. 

The LLE, j(N,u), is defined [H, M, HE US HI as 
the finite- TV version of Eq. ([9]). For the sake of com- 
pleteness, before going further we exhibit the evolution 
of the statistics of 7 (AT, lu) (hereafter referred to as 7jv for 
short) against increasing N. We consider the TM given 
in Eq. ®, with p = 1/2 and lu = 0.015 (fixed), for which 
case existing numerical results [Lsl. I20L l2l| give A w 57 
lattice spacings, i.e., 7 w 0.0175. 

In Fig. [H one sees that the probability distribution 
function (PDF), P(jn), of 7^ takes on an approximately 
Gaussian shape only for N > A . Fig. [2] (a) provides 
a quantitative check on the increasing narrowness of 
the PDF for N ^> A, and on the convergence process 
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Figure 1: (Color online) Normalized histograms of occurrence 
of LLE 7at for p = 1/2 (spin glass), for several values of N. 
For each histogram, N s = 10 samples were taken. Here, 
u = 0.015. 



0.03 




N -l/2 

Figure 2: (Color online) Estimates from the statistics of LLE 
7jv data, for N = 2 fe x 100, fc = 1, ••• ,9, against N~ 1/2 . 
For each AT, iV s = 10 5 samples were taken. Here, p = 1/2, 
io — 0.015. (a) Mean values {jn} and rms deviations An- 
Arrow on vertical axis indicates estimate from a single sample, 
with N — 10 7 . Full lines are, respectively, quadratic ((jn)) 
and linear (Ajv) fits of data, the latter taking only N > 800 
into account, (b) Skewness Sn- Error bars correspond to 
incompletely-sampled Gaussian distributions with N s = 10 5 
(see text). 



(in) — ► 7- One can infer the accuracy of the single- 
sample estimate corresponding to N — 10 7 (denoted by 
an arrow pointing to the vertical axis in Fig. [2] (a)), by 
extrapolating the TV-dependence of the width A^v of the 
PDF against N < 51200 (shown in the main diagram). 
The result is 7 = 0.01752(4), hence the associated error 
bar would be invisible on the scale of the Figure. 

The simplest quantitative indicator of whether the 
PDF actually turns Gaussian with increasing N [in the 
process of becoming a 5— function, characteristic of the 
self-averaging property just demonstrated] is its skew- 
ness [3l|, Sn = (((tjv — (jn))/An) 3 ). One must bear 
in mind that, because the number of samples N s is fi- 
nite (incomplete sampling), Sn itself will have a distri- 
bution. While this is true also for (7^) and An, the 
effects on Sn are much more prominent since it is the ra- 
tio of two quantities which, in the present case, both van- 
ish as N, N s — > 00. Indeed, we have checked that, for the 
data shown in Fig. [2] (a), corresponding to N s = 10 5 , esti- 
mates of (7jv) and An from distinct sequences of pseudo- 
random numbers usually differ only by a few parts in 
10 4 . On the other hand, it is known [31| that, for an in- 
completely sampled Gaussian distribution, the wi dth of 
its own skewness distribution is approximately \f&/N a 
(~ 0.0078 here). Since our own large- N PDFs are, to 
zeroth-order, Gaussian, we shall use this value as a lower 
bound to the uncertainty of our estimates for Sn- The 
results are displayed in Fig. [2] (b). Given the trend ex- 
hibited by the calculated Sn for the three largest values 
of N used (especially in contrast with the smaller- N re- 
gion), it appears safe to conclude that the skewness is 
in fact approaching zero, within the pertinent error bars. 
Thus, the overall evidence is compatible with a limiting 
Gaussian form for the PDF of the LLEs as N — > 00. 



B. Local effects and the Fokker-Planck equation 

We now turn to the opposite limit, which is our main 
concern here. Instead of considering the aggregate effect 
of contributions to 7jv, we shall analyse the PDFs of such 
contributions separately. 

The case N — 1 of Eq. JJJ]) is often denoted as "instan- 
taneous" LLE [28]. Here, we shall use the term in the 
following way. The instantaneous LLE at i = M, to be 
denoted by s(M), is the local contribution to the LLE, 
given at i = M, of the multiplying process denoted in 
Eq. 1J9J. One then has: 

1 N 

M=l 

For a one-dimensional mapping, the instantaneous LLE, 
as defined above, would be simply the local stretch fac- 
tor [2^]. In the present case, where even in Id the TM 
is 2 x 2, it is known that the eigenvectors (Lyapunov 
basis) of a product of random matrices are only local 
properties, as opposed to the corresponding eigenvalues 
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(i.e., the Lyapunov spectrum) which are global ones [23J]. 
Thus, one has a rotation of the local Lyapunov basis as 
the TM is repeatedly iterated. As the TM is symplec- 
tic, in Id this means that the (two) local eigenvalues are 
inverse of each other, so s(M) can be extracted by suit- 
able analysis of the growth factor associated to a given 
member of the Lyapunov basis. 

The connection with the F-P equation proceeds as fol- 
lows. One gets, by taking the continuum limit of Eq. (7|): 



u((£)u{£) 



d 2 u{t) 



(11) 



where £ stands for position along the axis. One sees that 
u(€) is the result of accumulated contributions from the 
instantaneous Lyapunov exponent s(£'), < £' < I: 



u(£) =exp / s(£')d£' 



(12) 



Therefore, the right-hand side of Eq. fTTj) turns into: 



d 2 u(£) 



* 2 (t) + ^f) <«"- (13) 



Eq. (fTTj) then yields the following equation for the 
stochastic variable s{£): 



(14) 



In the general case £ is a binary-distributed variable with 
mean M.\ = 2p — 1 and variance M.2 = 4p(l — p). The 
F-P equation corresponding to the evolution of the prob- 
ability distribution of s is^with £ naturally assuming the 
role of time- like variable llfl : 



dP(sJ) _ d 



Of 



ds 



((s 2 +uMi)P(s,e)) + 
-\^ 2 {M^P{s,£)) . 



(15) 



C. The unmagnetized HMSG 



In this subsection we take the special case of an un- 
magnetized HMSG (p =1/2). Eq. (HU then becomes: 



dPjs, 

ae 



^Pi S J)) + \^ 2 ^ 2p M) ■ (16) 



Assuming stationarity, P(s,£) — > P(s) (which corre- 
sponds to the regime £ > A), the equation to solve is: 

|(. s F(»)) + i^(^F(»))=0. (17) 



Thus, one gets: 

P(s) = C exp (-2s 3 /3w 2 ) 



exp (2y 3 /3w 2 ) dy . 



(18) 
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Figure 3: (Color online) Normalized histograms of occurrence 
of instantaneous LLE s(N) for several values of N. For each 
histogram, iV s = 10 6 samples were taken. Here, p = 1/2, 
lu = 0.015. 



It is immediate to see that the 3/2 power, referred to 
in connection with Eq. J6]), arises in the scaling variable. 
Furthermore, it is clear from the derivation of the present 
results that they apply only in the scaling regime, i.e., 
u -> 0, small |s|, £ > A. 

With dy dx exp(x 3 - y 3 ) = C w 4.8, the 



normalization of Eq. lfT8|l implies 



C = C(u) 



1 



(3w 2 /2) 1 /3c 



(19) 



The numerical verification of the ideas just presented 
begins by checking the behavior of the PDFs for s(N), 
against varying AT. Selected data are shown in Fig. [31 
One can see that, for the ranges of N considered in the 
Figure, the overall shape of the PDFs remains roughly 
constant, though a slight narrowing and shifting of the 
central peak take place as N increases. This is in contrast 
with the LLEs exhibited in Fig. [H for which variation of 
N in the same range was accompanied by a pronounced 
change in shape and width of the respective PDFs. For 
1 < N < 13 [not shown in the Figure], the PDF for 
s(N) starts nearly flat at N = 1, and gradually devel- 
ops both the central peak, and the small "wings" at the 
ends (whose relevance in Fig.[3jis overemphasized by the 
logarithmic scale on the vertical axis). The very good 
superposition of the N = 50 and 200 data in the Fig- 
ure indicates the predicted trend towards a fixed, TV- 
independent form, which corresponds to stationarity of 
the F-P equation [see Eq. (fT7|)J. This has been confirmed 
by examination of numerical PDFs for 100 < TV < 600. 
Recalling that the localization length for w = 0.015, as 
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Figure 4: (Color online) Squares: normalized histogram of 
occurrence of instantaneous LLE s, from N s — 10 6 samples, 
for N = 200, p = 1/2, to = 0.015. Continuous line is the 
analytical form, Eq. (JHJ, normalized and with appropriate 
scaling of variables [see also Eq. I|19p ], 



Figure 5: (Color online) Scaling plots of normalized his- 
tograms of occurrence of instantaneous LLE , as suggested 
by Eq. (| 1 8 j) . for assorted energies. In all cases, p = 1/2, 
iV s = 10 6 samples, N = 200 (so the condition N > A(«) is 
always obeyed). 



is the case of Fig. 02 is A w 57 lattice spacings, we see 
that the condition for an iV-independent PDF is indeed 
N > A. 

We now test whether the form given in Eq. lfT8|) de- 
scribes our numerical data. Figure [4] shows that the 
agreement is excellent, in the region —0.5 < s < 0.5, 
which is where scaling is expected to hold. Note that 
there are no adjustable parameters in the fit, all scale 
factors being provided by normalization and scaling con- 
siderations. 

In Fig. [5] we test for consistency of scaling among nu- 
merical data for different energies w, as suggested by 
Eq. l[T8|) . Again, the agreement with predictions is excel- 
lent, provided that the conditions upon which Eq. lfl8|) 
was derived are obeyed. 



D. The magnetized HMSG (p §) 

For the magnetized case p / 1/2 we return to the gen- 
eral form of the F-P equation, Eq. lfl~5|). whose stationary 
solution is given by: 

P(«) = C'exp(-/( fl )) f exp(f(y))dy, (20) 

J — DO 

where 

f( v )=f(y,cj,p) = -± r (^ + 2LuM iy ) . (21) 



In analogy with Sec. lIII C| we first check whether Eq. f20|) 
gives a faithful description of simulational data, for given 
p and to. In Figure 03, calculated data for p = 0.8, 
lu = 0.015 are shown, together with the analytical re- 
sult from Eq. f20|) for the same values of p and lu. Apart 
from an overall normalization factor which affects only 
the vertical scale, there are no adjustable parameters. 
Again, the agreement is very good. 

In order to emphasize the degree to which the shape 
and range of the PDFs are affected by variations in the 
physical parameters, Figure [6] also shows the PDF for 
the unmagnetized case, at the same energy w = 0.015. 
As p increases from 1/2, the general trend is towards 
reduction of the skewness (which is w —0.32 at p = 1/2, 
—0.02 at p = 0.80, and [in absolute value] < 0.01 at 
p = 0.95). The peak of the distribution, which is about 
one standard deviation away from the origin at p = 1/2, 
approaches s — as p approaches unity. This process is 
accompanied by a broadening of the PDF. 

We now proceed to framing the preceding observations 
within a crossover description. One sees that, because 
f(y) in Eq. l(2T|) is a polynomial in y, for p ^ 1/2 it 
is not possible to develop scaling arguments (over the 
whole range of the variable s) similar to those invoked 
for p — 1/2, and graphically depicted in Fig. [H Nev- 
ertheless, analysis of Eq. J21|) shows that the crossover 
away from unbiased (p = 1/2) behavior is governed 
by the dimensionless ratio z = (lj Mi) 1 / 2 /^ 2 A4 2 ) 1/3 
(w [(p - l/2)/w 1 / 3 ] 1 / 2 away from p -> 1). This way, s 
scales with lj 2 / 3 for z <§; 1, and with u/(p — 1/2) for 
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Figure 6: (Color online) Hexagons: normalized histogram of 
occurrence of instantaneous LLE s, from iV s = 10 6 samples, 
for to = 0.015, p = 0.80 (here we used N = 600, as the local- 
ization length is A » 498). Continuous line is the analytical 
form, Eq. I|20p . normalized. The dashed line is the analytical 
result for p = 1/2, from Eq. lfl8l) . 

z> 1. Analysis of the asymptotic behavior of P(s), as 
given by Eqs. 1|20|1 and ([2T) . shows that, with z defined 
as above and v = suj~ 2 / 3 , one has for the tails of the 
distribution: 

R(v) ~ -^-^ \v\ » 1 , (22) 

where i?(v) = w 2 / 3 P(s), and (for z ^> 1) the normaliza- 
tion constant C oc z. 

Numerical checks of Eq. f22|) against simulational data 
must be carried out for suitable ranges of p, u), s (equiv- 
alent^, z, v), such that (i) asymptotic behavior has al- 
ready set in (i.e. large \s\, \v\), and (ii) one is still within 
the scaling regime (which implies small \s\). Thus one 
expects Eq. (|22| to hold only within windows of varying 
width, which may or may not be easy to identify against 
statistical noise. 

In Figure [7] we show that it is indeed possible to find 
intervals of v along which a plot of (v 2 + z 2 ) R(v) is 
rather flat, as suggested by Eq. lj22j) . While the evidence 
is somewhat smeared for the upper set of data (corre- 
sponding to p = 0.95, u> — 0.015, i.e. z\ « 3.32), it is 
clear for the lower one (for which p = 0.80, w = 0.0075, 
i.e. 22 ~ 2.03). Furthermore, the ratio of the averages 
of (v 2 + z 2 ) R(v) (each denoted by a horizontal line in 
the Figure) along the respective flat sections [these lat- 
ter defined with an inevitable degree of arbitrariness] is 
« 1.4(1). This compares reasonably well with the pre- 
dicted value Z1/Z2 — 1.64, see Eq. ([22| and the comments 



Figure 7: (Color online) Symbols are data from normalized 
histograms of occurrence of instantaneous LLE s, from N s = 
10 6 samples, plotted as suggested by Eq. (1221) . See text for 
definitions of z, v, R(v). Hexagons: p = 0.8, lo = 0.0075 (z w 
2.03); squares: p = 0.95, u = 0.015 (z w 3.32). Horizontal 
lines mark central estimates of averages of (v 2 + z 2 )R(v), 
respectively along —9 < v < —3 (lower) and —8 < v < —4 
(upper). 



immediately below it. Given that the latter result is ex- 
pected to hold for z ^> 1, one may infer that the small 
discrepancy found is due to Z\, Z2 not being large enough. 

Finally we remark that the pure-ferromagnet limit is 
somewhat subtle, because in order to reach the station- 
ary regime of the F-P equation where the iV-independent 
PDF forms hold, the condition N > X must be obeyed. 
However, as p — > 1, and for the low energies where scaling 
is valid, A diverges as one approaches the pure-system 
magnon band (located at < ui < 2, in the current 
units). For example, at p — 0.95 and u> = 0.015, one gets 
A w 2500. 



E. Non-stationary regime for p 7^ | 

Still for the magnetized case p / 1/2, one can 
learn more by considering selected aspects of the non- 
stationary regime. Going back to Eq. (fl5| . one sees that, 
in the limit p— 1/2 3> oj 1 / 3 , i.e., ?»1, the diffusive term 
becomes small and this regime is associated with just 
streaming in leading order. The non-stationary equation 
to be solved is thus: 

^^M^^+^OPM) ■ (23) 
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Figure 8: (Color online) Central sections of normalized his- 
tograms of occurrence of instantaneous LLE s, from N B = 10 7 
samples, for uj = 10 -6 , p = 0.60 and N as indicated. 



Figure 9: (Color online) Scaling plot of data displayed in 
Figure [H as suggested by Eq. I|25p. 



With P = [to (p - 1/2)] 1 /2 ) t =l/3,y = s/0, and defining 
Q(t, y) such that Q(t,y)dy = P(s,£)ds, Eq. $F5}) turns 
into: 



-dt=d-y( y 



whose general solution is: 



1 



Q(i,y) = -^/(t + tan^y) 
1 + y 



(24) 



(25) 



In order to have the condition p-l/2> uj 1 ^ 3 fulfilled to 
a good extent, so that the scaling behavior predicted by 
Eq. lf25j) could be unequivocally demonstrated, we used 
p = 0.6, uj = 10 -6 . For such a choice, the localization 
length is A > 4 x 10 . This allowed us to take values of 
N equal to several hundred lattice spacings, which are 
both much larger than unity (so discrete-lattice effects 
are negligible) , and still much shorter than A (thus guar- 
anteeing non-stationarity by a broad margin). Figure [8] 
shows the region close to the central peaks of the PDFs 
(where scaling is expected to hold), for N — 150, 250, 
500 and 750. In Fig. [9] the same data are shown in a 
scaling plot, as suggested by Eq. f25|) . Note that the 
number of samples N s is one order of magnitude larger 
than, e.g., that used in earlier sections of this work. This 
was necessary, in view of the relatively wide scatter of the 
PDFs: typically, they displayed almost flat tails running 
out to \s\ > 0.07, so the relevant data as far as scaling is 
concerned were a small subset of the total gathered (com- 
pare the horizontal scale in Fig. [8]). Even so, one can see, 
close to the bottom of the scaling curve in Fig. [9j that a 



non-negligible degree of fluctuation-induced spread still 
remains. Nevertheless, overall agreement with the scaling 
predictions of Eq. (|25|) is remarkable. 



IV. DISCUSSION AND CONCLUSIONS 

We have investigated scaling properties of the probabil- 
ity distribution functions (PDF) of Lyapunov exponents 
for the one-dimensional Heisenberg-Mattis spin glass. 

In Section UlI Al we showed that, for a given energy uj 
(small enough, such that the localization length A(w) is 
sufficiently large for scaling concepts to apply) , the PDF 
of the local Lyapunov exponent (LLE) j(N,ui), takes on 
a shape which, for N > A(w), is increasingly close to 
a Gaussian. In the limit TV — > oo, the PDFs turn into 
5 functions, in conformity with the central limit theo- 
rem. Thus, such aggregate effects reflect only general sta- 
tistical properties of fluctuations in systems with many 
(almost-) independent degrees of freedom. The connec- 
tion with the underlying physical problem is traced ex- 
clusively through the dependence of the numerical value 
of the (asymptotic) Lyapunov exponent with energy (and 
ferromagnetic bond concentration; though this latter as- 
pect was not exploited here, it has been investigated be- 
fore [is, mm). 

On the other hand, the PDFs for the instantaneous 
Lyapunov exponents, as defined via Eq. (fTOf . display 
the following properties of interest, exhibited in Sec- 
tions ITTTBl mfcl and ITTTdI (i) for N > A(w) they 
approach a fixed, non-trivial shape with nonvanishing 
width; (ii) such specific shape can be predicted through a 
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connection with the stationary state of a Fokker-Planck 
(F-P) equation; and (iii) the shape is a "fingerprint" 
shared by all the systems in the universality class in which 
the statistical properties of the system lie . Furthermore, 
the F-P equation providing the link can be set up directly 
from, and closely reflects, the physical features of the sys- 
tem under investigation [ see especially Eqs. lfTTj) - (fT6"l) ]. 

The applicability of a description via an F-P equa- 
tion goes beyond the stationary state, as shown in Sec- 
tion MI El There, it is shown that the scaling of non- 
stationary PDFs in a specific regime (in which diffusion 
effects are expected to be negligible) closely follows pre- 
dictions drawn from the corresponding form of the F-P 
equation. 

In the present work we have demonstrated that an F- 
P approach can be successively applied to both station- 
ary and non-stationary properties of the PDF of instan- 
taneous Lyapunov exponents for Heisenberg-Mattis spin 
glasses. 

As a final remark, we believe that the scaling proper- 
ties of the PDF of Lyapunov exponents described here 
are identical to those pertinent to the zero-temperature 
random-bond (±1) classical Heisenberg chain. This 
comes via the identification of the scaling properties of 
the low-energy excitations of the HSMG chain with those 



of the classical Heisenberg model in Id flR ana a 
plus the general applicability of the F-P equation in the 
continuum 

It is expected that treatments along these lines can be 
devised, with similar degree of success, for other physical 
systems in whose description Lyapunov exponents play a 
prominent role. 
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